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Abstract 



Magnetization and Cu NMR-NQR relaxation measurements are used to 
study the superconducting fluctuations in YBa2Cu306+a' (YBCO) oriented 
powders. In optimally doped YBCO the fluctuating negative magnetization 
Mfi(H,T) is rather well described by an anisotropic Ginzburg-Landau (GL) 
functional and the curves Mji/^/H cross at Tc. In underdoped YBCO, in- 
stead, over a wide temperature range an anomalous diamagnetism is observed, 
stronger than in the optimally doped compound by about an order of mag- 
nitude. The field and temperature dependences of Mji cannot be described 
either by an anisotropic GL functional or on the basis of scaling arguments. 
The anomalous diamagnetism is more pronounced in samples with a defined 
order in the Cu(l)0 chains. The ^'^Cu(2) relaxation rate shows little, if any, 
field dependence in the vicinity of the transition temperature Tc{H = 0). It 
is argued how the results in the underdoped compounds can be accounted for 
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by the presence of charge inhomogeneities, favoured by chains ordering. 
PACS numbers: 76.60.Es, 74.72.Bk, 74.20.Mn 
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I. INTRODUCTION 



The superconducting transition in conventional superconductors is rather well described 
by mean field theories, essentially because in the coherence volume a large number of 
pairs is present. On the contrary, in cuprate superconductors the small coherence length 
C,, the reduced carrier density, the marked anisotropy and the high transition temperature 
Tc strongly enhance the superconducting fluctuations (SF). In a wide temperature range, 
which can extend up to 10 or 15 K above Tc, a variety of phenomena related to can 
be detected by several experiments such as specific heatS, thermal expansion!, penetration 
depthi, conductivityB and magnetization! measurements. Also NMR-NQR relaxation has 
been used to detect Recently the field dependence of the NMR relaxation rate W 

and of the Knight shift have been studied, by varying the magnetic field from 2 up to 24 
Tesla0§. Mitrovic et al.tll attributed the field dependence of W{H) to the corrections in 
the density of states (DOS) contribution to the spin-lattice relaxation. On the other hand, 
Gorny et al.El have observed W{H) field independent, up to 14 Tesla, from 150 K to 90 K, 
in a sample showing a decrease in 1/TiT starting around 110 K. An analysis of the role of 
SF in NMR experiments and of the field dependence of W has been recently carried out by 
Eschrig et al.0, by extending the analitical approach of Randeria and Varlamovlii to include 
short wave-length and dynamical fiuctuations. 

Of particular interest is the role of SF in underdoped high temperature superconductors. 
In fact, in these compounds one has the peculiar phenomena of the opening, at T* ^ Tc, of 
a spin gap (as evidenced by NMR and neutron scattering) and of a pseudo-gap (as indicated 
by transport and ARPES experiment These gaps, which are possibly connectedS, 

have been variously related to SF of particular character&iH. 

In this paper magnetization and ^^Cu{2) NMR-NQR relaxation measurements in under- 
doped YBCO are reported and compared with the ones in the optimally doped compound. 
The magnetization data yield information on the fiuctuating diamagnetism (FD), while the 
nuclear relaxation rates W convey insights on the effects of SF on the k-integrated spin 
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susceptibility in the zero frequency limit. 

The paper is organized as follows. In Sect. II some basic equations describing SF and 
FD are recalled. After some experimental details (Sect. II), in Sect. Ill the experimental 
findings and their analysis are presented, with emphasis on the effect of the field on W and 
on the behavior of the diamagnetic magnetization in chain-ordered underdoped YBCO. The 
main results and conclusions are summarized in Sect. IV. 



II. THEORETICAL FRAMEWORK AND BASIC EQUATIONS 

Because of SF the number of Cooper pairs per unit volume, which is given by the average 



value of the square of the modulus of the order parameter y< {ipl"^ >, is different from 
zero above T^. In the time-dependent Ginzburg-Landau (GL) descriptionlli the collective 
amplitudes and the correspondent decay times of SF are given by 

< l^-^l >= TTe^ = YTW' 

where tgl = (vr^SfesTje-^ and ^(T) = ^{0)e-^^^ , with e = (T - Tc)/Tc. < l^-kP > plays 
the role of the Fourier components of the average number of pairs per unit volume, while Tk 
is the correspondent relaxation time. 

Fluctuating pairs can give rise to a diamagnetic magnetization above Tc. The diamag- 
netic magnetization Mfi in general is not linear in the field H. The curves Mfi vs. H and 
Mfi vs. T can be tentatively analized by generalizing the Lawrence-Doniach model, using 
the free energy functionaliiS 



4m|| nc 



(2) 



where the last term takes into account the tunneling coupling between adjacent layers. From 
Eq. 2, deriving the free energy in the presence of the field and by means of a numerical 
derivation {Mfi = —dF/dH), one can obtain the fluctuating magnetization. 



According to scaling arg umentHhi for moderate anisotropy (quasi-3D case) one expects 
that the magnetization curves, at constant field, cross at Tc{H = 0) when the magnetization 
is scaled by \/H. The amplitude Mji at T = Tc departs from the Prange's result@ by 
a factor 3 — 7, corresponding to the anisotropy ratio Accordinglyii, the data at 

different fields collapse onto an universal curve when Mfi/T is reported as a function of 
e/H^, for a critical exponent for the coherence length corresponding to ~ 0.66. 

For strong anisotropics, namely quasi-2D systems, the Mfi curves at constant field cross 
each other at Tc(0). Mfi{Tc) is larger than the value obtained from the Gaussian approxi- 
mation by a factori around 2. Collapse of the data onto a universal curve occurs for u well 
above 0.66, usually in the range 1.2 — 1.4. 

The contributions to the relaxation rates W due to SF can be derived within a Fermi 
liquid scenario, withouth specifying the nature of the interactionsi. The direct and most 
singular contribution, equivalent to the Aslamazov-Larkin paraconductivity, is not effective 
as nuclear relaxation mechanism. The positive Maki- Thompson (MT) contribution Wmt 
results from a purely quantum process, involving pairing of the electron with itself at a 
previous stage of motion, along intersecting trajectories. A negative SF contribution Wdqs 
comes from the density of states reduction when electrons are subtracted to create pairs. 
Approximate expressions for Wdos and Wmt can be derived by resorting to simple physical 
arguments as follows. The relaxation rates can be approximated in the form 

2W ^ ^{A^)ksTY. ( x"sak,-^h \ ^ ^(^A^^ksTxsMO,0)T.r^, (3) 

where < > is an average form factor for Cu nuclei and = Jg(0) an effective spectral 
density0@. SF modify the static spin susceptibility Xspm(0, 0) and the effective correla- 
tion time J^k'^k From Eq.l the in-plane density of pairs is ric = Z^fcd^gP) = 
2ne{kBTc/ Ep) In ^ and therefore the SF imply 

(4) 
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For the dynamical part in Eq. 3, the DOS contribution can be obtained by averaging over 
the BZ the collective correlation time in Eq. 1: 



{tdos) = 



(5) 



in the dirty limit, differing from the exact calculationlli only by a small numerical factor, r 
is the electron collision time. The MT contribution has to be evaluated in the framework 
of diagrammatic theoriesil. The final result corresponds to the 2D-average over the BZ of a 
decay rate of diffusive character = Dk/^ ( with D = carrier diffusion constant), 

phenomenologically accounting for phase breaking processes by adding in the decay rate a 
frequency r^"^: 



{rMT) = T.{iDk' + r^Ml + ek')}-' 

k 



1 



In- 



(6) 
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where 7,^ = ^'^{Q)/t^D = uTi / 9>K bTcT^ is a dimensionless pair breaking parameter. 

By indicating with the relaxation rate in the absence of SF and by neglecting the 
correction to Xspin(0,0) (Eq. 4), from Eqs. 3, 5 and 6 the relaxation rate in 2D systems 
turns out 



W 



SF 
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^EpT (e - 7,^) 70 EpT e 
To extend this equation to a layered system, when dimensionality crossover (2D^ 3D) oc- 



curs, one has to substituteB /n(i) by 2/n 1l{\ft + Ve + r) and /n(;^) by 2/n 



_(VhV£±D_ 



with r = 2^c^(0)/c?^ anisotropy parameter [d is the interlayer distance). It is noted that the 
MT contribution (first term in Eq. 7) is present only for s-wave orbital pairing, while it is 
averaged to almost zero for d-symmetry. 



III. EXPERIMENTAL DETAILS 

The measurements have been carried out in oriented powders of optimally doped YBCO 
in one chain-disordered underdoped and in two chain-ordered underdoped YBCO sample 
The oxygen content in the underdoped compounds was close to 6.66, with sligth differences 
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in Tc- Electron diffraction microscopy evidenced the expected tripling of the a-axis, while 
resistivity measurements show a sharp transition with zero resistivity at T= 62 K and 
occurrence of paraconductivity below about 75 Table I collects the main properties of 
the samples, as obtained from a combination of measurements. 

The ^^Cu relaxation rates 2W = Ti~^ have been measured by standard pulse techniques. 
In NQR the recovery towards the equilibrium after the saturation of the (±| ±|) line is 
well described by an exponential law, directly yielding GW . For the NMR relaxation a good 
alignment of the c axis of the grains along the magnetic field is crucial to extract reliable 
values of W (H). The NMR satellite line, corresponding to the (| |) transition can be 
used to adjust the alignment and to monitor the spread in the orientation of the c-axis, the 
resonance frequency being shifted at the first order by the term cQV^z^'icos'^O — l)/h, due 
to the quadrupole interaction ( 6 angle of the c-axis with the field). From the width of the 
spectrum (Fig. la) the spread in the orientation of the c-axis appears within about two 
degrees. 

The recovery law for the NMR satellite transition is 

y{t) = 0.5exp{-12Wt) + 0Aexp{-6Wt) + 0.1exp{-2Wt) (8) 

This law has been checked to fit well the experimental data (Fig. lb), proving the lack of 
background contamination due to other resonance lines. 

The magnetization has been measured by means of a Metronique Ingegnerie SQUID 
magnetometer, on decreasing the temperature at constant field and, at selected tempera- 
tures, by varying H. The paramagnetic contribution to M was obtained from M vs if at T 
> 110 K, where practically no fluctuating magnetization is present. Then Mfi was derived 
by subtraction, singling out a small contribution from paramagnetic impurities. The param- 
agnetic susceptibility turns out little temperature dependent around and this dependence 
will be neglected in discussing the much stronger diamagnetic term. 

In optimally doped YBCO the measurements have been carried out at the purpose 
to confirm the results recently obtained in single crystals by other author sUSiiii. The 
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isothermal magnetization curves are satisfactorily described on the basis of the anisotropic 
GL functional ( Eq. 2), for e > 4 ■ 10~^. Only close to Tc and for H < 1 Tesla an observable 
departure is detected, indicating crossover to a region of non-Gaussian SF, in agreement 
with recent thermal expansion measurementsi. The Mji vs. T curves cross at Tc(0) when 
Mfi is scaled by y/H, as expected for moderate anisotropy0. 

IV. RESULTS AND ANALYSIS 

Let us first comment the data in optimally doped YBCO (Fig. 2a). The comparison 
between the NMR and NQR W^s has been already discussed in previous papersi'0. Here we 
only add a few comments motivated by more recents worksO'ii involving the remarkable 
aspect of the field dependence. 

The NQR W can be compared to Eq. 7, by using the value^l r = 10~^^s, = 2.10~^'^s 
and a dephasing time parameter 7^ = 0.2. A quantitative fitting is inhibited by the fact that 
the background contribution to W does not follow the Korringa law, in view of correlation 
effects among carriers^!. A firm deduction, however, is that an anisotropy parameter r 7^ 0, 
and namely a crossover to a 3D regime, is required to avoid unrealistic values (2D line in 
Figure 2a), as indicated also by magnetization measurements (see later on). 

Mitrovic et al.0 have discussed the field dependence of Ti in terms of quenching of the 
DOS term only, by resorting to the theory of Eschrig et al.0. The field dependence of the 
SF contribution to the relaxation rate is a delicate issue, because of the nontrivial interplay 
of many parameters which include reduced temperature e, reduced field P = j^, anisotropy 
parameter r, elastic and phase breaking times r and r^. Eschrig et al.0 have extended 
the analitical approach^ by taking into account arbitrary values of [KbTt short-wave 
fluctuations and dynamical fluctuations. The price of this generalization is the restriction to 
a 2D spectrum of SF fluctuations. In view of our experimental flndings the analysis of the 
fleld effect in a purely 2D framework is questionable. On the other hand, a re-examination of 
the field effect for a layered system is now possible by applying the method of the DOS term 
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regularization devised by Budzin and Varlamovil, which has indicated how the divergences 
can be treated. Since the dynamical fluctuations could be relevant only for fields of the 
order of ifc2, up to 20 — 30 Tesla the fluctuations can be safely treated in the nearly static 
limit. The fleld dependence for weak flelds {(3 <^ e) turns outS 

Wi, « e) - WiO, .) ^ W^-<»'^ (.(Tr) - ^) ^l±l|^,^ (9) 

where 

ri4C(3)/7r3, Tr«l 



k{Tt) 



7C(3) 1 

TT 47rTr /2)-i/'(l /2+1 /4ttTt)] +ip'{l/2) 



ATt , K Tr < 1 



(here h = ks = 1). 

So both effects of increase and of decrease of W on increasing the magnetic fleld are 
possible, depending on the mean free path in the speciflc sample. If Tr ^ 0.1 the main 
correction to W is due to the MT term and one should observe W decreasing with increasing 
field, while if Tr ^0.1 the DOS correction becomes dominant and W is expected to increase. 

As it appears from Fig. 2 the effect of the fleld for T > Tc{H = 0) is small, if any. It 
is noted that if Eq. 7 is applied to the underdoped compounds the reduction in the Fermi 
energy Ep and the increase in the anisotropy (i.e. decrease of r) would imply a sizeable 
increase in W^^. This enhancement in the underdoped compounds does not occur (see 
Fig. 2b and 2c for the chain-disordered and chain-ordered YBCO's, respectively), unless 
the decrease of W over a wide temperature range, usually related to the spin-gap opening0, 
should be attributed to a fleld-independent DOS term. On the other hand, the conventional 
SF of GL character should occur only close to T^, where, on the contrary, Xspm(0, 0) is little 
temperature dependent. Eq. 9, in principle, does predict a fleld dependence. However, 
if typical values Ep ^ 3000 K, r ~ 10"^'', {k{Tt) — (7r/87<^) ~ 1 and r ~ 0.1 are used, 
for e ^ 3 X 10"^ one has W{(3) - W{Q) ^ O.ISPF^^^/^^, hardly to evidence for /3 < 1. A 
comprehensive discussion of the fleld dependence of the SF contribution to the relaxation 
rate is given elsewhere^. Here we only remark that the results in strong flelds {(3 ^ 0.2) 
from Mitrovic et al.0 can hardly be justifled on the basis of Eq. 9. 



As regards the magnetization measurements in optimally doped YBCO, as already men- 
tioned, our data indicate a 3D regime crossing from Gaussian to critical fluctuations close 
to Tc- The value of 7a„m3(oo)0 at the crossing point of the curves Mfi/y/H vs. T turns out 
around —1.5, corresponding to an anisotropy factor = 4.5 . A collapse onto a common 
function is obtained when Mfi/\fHT is plotted as a function of e/H^''^^'^ . The collapse fails 
in magnetic fields less than about 0.3 Tesla. A small-field departure from the universal 
function has been already observed in underdoped compoundsil, but no mention about it 
is found in the literature for optimally doped YBCO. 

In the underdoped compounds a marked enhancement of the fluctuating magnetization 
is observed (Figs. 3 and 4). In Fig. 4 the susceptibilty, defined as ^ ioi H = 0.02 Tesla, is 
reported. The inset is a blow up of the results around the temperature where reversing of the 
sign of the magnetization occurs. The susceptibilities for chain-ordered and chain-disordered 
underdoped samples are compared with the one in optimally doped YBCO in Fig. 5. 

In underdoped YBCO the magnetization curves dramatically depart from the ones ex- 
pected on the basis of Eq. 2 and numerical derivative. No crossing point is observed in the 
curves Mji vs T and no collapse is found on anisotropic 3D or 2D-like curves for Mfi/\/H 
or Mfi vs tjH^I'^^ (Fig. 6). These anomalies are more manifested in the chain ordered 
compound, as it can be realized from the peculiar shape of the isothermal magnetization 
vs. H (Fig. 7). One could remark that a fraction of a few percent of non-percolating local 
superconducting regions can account for the absolute value of the diamagnetic susceptibility 
and for the shape of the magnetization curves. A trivial chemical inhomogeneity could be 
suspected. The fact that the anomalous diamagnetism has been observed in three samples 
differently prepared and since it is strongly enhanced in the chain-ordered ones, supports 
the intrinsic origin of the phenomenon. Recent theorieJ^lS have discussed how a distribu- 
tion of local superconducting temperatures, related to charge inhomogeneities, can cause 
anomalous diamagnetic effects. 

Mesoscopic charge inhomogeneities have been predicted on the basis of various theoreti- 
cal approacheJnilii. In underdoped compounds one could expect the occurrence of preformed 
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Cooper pairs causing "local" superconductivity lacking of long-range phase coherence. An 
inhomogeneous distribution of carriers, on a mesoscopic scale, is supported by ab plane opti- 
cal conductivity^ and it is the basic ingredient of the stripes modefl. Since the anomalous 
diamagnetism is enhanced in chain-ordered underdoped YBCO one is tempted to relate it 
to the presence of stripes. The insurgence of phase coherence among adjacent charge-rich 
regions can qualitatively be expected to yield strong screening above the bulk superconduct- 
ing temperatureiiiiS. Therefore, one could phenomenologically describe the role of SF in 
underdoped compounds, as follows. At T ^ Tc, just below T*, fluctuating pairs are formed, 
without long-range phase coherence. Below Tsf, about 20 — 30 K above Tc, phase coherence 
among adjacent charge-rich dropsS start to develop, yielding the formation of supercon- 
ducting loops with strong screening and causing the onset of the anomalous diamagnetism. 
Close to Tc, SF of GL character occur and long-range phase coherence sets in. 



V. SUMMARIZING REMARKS AND CONCLUSIONS 

By combining magnetization and NMR-NQR relaxation measurements, an attempt has 
been done to clarify the role of superconducting fluctuations and of fluctuating diamag- 
netism in underdoped YBCO, vis-a-vis the optimally doped compound. In the latter case 
the fluctuations above Tc are rather well described by an anisotropic Ginzburg-Landau (GL) 
functional and by scaling arguments for slightly anisotropic systems. A breakdown of the 
Gaussian approximation for small magnetic fields has been observed close to Tc. The ®^Cu 
relaxation rates W around Tc show little field dependence, if any. One cannot rule out the 
presence of a MT contribution, and then of a small s— wave component in the spectrum of 
the fluctuations, which should be sample-dependent in view of the role of impurities in the 
pair-breaking mechanisms. In underdoped YBCO an anomalous diamagnetism is observed, 
on a large temperature range. The diamagnetic susceptibility at Tc is about an order of mag- 
nitude larger than the one in the optimally doped sample and the isothermal magnetization 
curves cannot be described by the anisotropic GL functional. Scahng arguments, such as 
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the search of an universal function in terms of e/H^/'^^ , appear inadequate to justify the 
experimental findings. The anomalies are more marked in the chain-ordered samples. Also 
in the underdoped YBCO ^^Cu W in NQR almost coincide with the NMR ones aX H = 9.4 
Tesla, in agreement with a theoretical estimate of the effect of the magnetic field. Another 
conclusion drawn from the field and temperature dependences of the ^^Cu relaxation rate is 
that the spin gap does not depend on the magnetic field and that the behavior of the spin 
susceptibility for T ^ Tc cannot be ascribed to SF of GL character. 

The fluctuating diamagnetism observed in underdoped YBCO can be phenomenologically 
accounted for by the presence of charge inhomogeneities at mesoscopic level. The anomalies 
in the magnetization curves can be attributed to " locally superconducting" non-percolating 
drops, present above the bulk Tc and favoured by chain ordering. 
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FIGURES 

FIG. 1. a) NMR spectrum of the ^^Cu(2) satellite line and the correspondent relaxation law 
(6) (the solid line is the best fit according to Eq. 8 in the text). 

FIG. 2. a) NQR (o) and NMR (□) {Ho = 6 Tesla || c) '^^Cu relaxation rates 2W in optimally 
doped YBCO. The solid and dashed lines are guides for the eye. The dotted line is the temperature 
behavior according to Eq. 7 in the text for r = (namely 2D spectrum of fluctuations) and 7^ = 0.2. 
NQR (empty simbols) and NMR (full symbols) relaxation rates 2W (in a field of 9.4 Tesla along 
the c-axis) in chain-disordered (b) and chain-ordered (c) underdoped YBCO around Tc. 

FIG. 3. Raw data for M/H {H = 0.1 Tesla) in chain-ordered YBCO, for field along the c-axis 
( □ ) and in the a-b plane ( A ). In the inset the results for Hc2 around Tc are shown (estimated 
experimental error < 2 %). The solid lines indicate the procedure to evaluate Tc{H) and rc(0). 

FIG. 4. Ratio M/H, in a field of 0.02 Tesla along the c-axis, in chain-ordered YBCO. The 
inset is the blow-up of the data for T > 85K. Similar results have been obtained in another sample 
prepared with the same procedure. 

FIG. 5. Comparison of the susceptibility (defined as M/H, for H = 0.02 Tesla along the 

c-axis) in optimally doped YBCO (squares) and in chain-disordered (circles) and chain-ordered 
(triangles) underdoped YBCO. 

FIG. 6. Evidence of the failure of anisotropic GL functional and of scaling arguments to 
describe the magnetization curves in underdoped YBCO. The solid lines in the upper part of the 
figure are the behaviours obtained from the numerical derivative of Eq. 2. The dashed lines are 
guides for the eye. 

FIG. 7. Magnetization curve in chain-ordered underdoped YBCO under high stabilization 
temperature {T = 65.67K ). a) raw data for the total magnetization, b) fluctuating magnetization 
Mfi as a function of the field. 
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